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We consider bipartite systems as versatile probes for the estimation of transformations acting locally on one 
of the subsystems. We investigate what resources are required for the probes to offer a guaranteed level of 
metrological performance, when the latter is averaged over specific sets of local transformations. We quantify 
such a performance via the average skew information, a convex quantity which we compute in closed form for 
bipartite states of arbitrary dimensions, and which is shown to be strongly dependent on the degree of local purity 
of the probes. Our analysis contrasts and complements the recent series of studies focused on the minimum, 
rather than the average, performance of bipartite probes in local estimation tasks, which was instead determined 
by quantum correlations other than entanglement. We provide explicit prescriptions to characterize the most 
reliable states maximizing the average skew information, and elucidate the role of state purity, separability and 
correlations in the classification of optimal probes. Our results can help in the identification of useful resources 
for sensing, estimation and discrimination applications when complete knowledge of the interaction mechanism 
realizing the local transformation is unavailable, and access to pure entangled probes is technologically limited. 
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I. INTRODUCTION 

Quantum metrology is one of the most promising branches 
of quantum technology and studies how to exploit the laws 
of quantum mechanics to improve the precision in the estima¬ 
tion or identification of some target parameter characterizing 
a quantum system of interest [1-5]. A typical estimation sce¬ 
nario involves three distinct phases [3]: (i) a probe system is 
initialized in an input state; (ii) the probe interacts with the 
system that encodes the parameter to be estimated; (iii) the 
output state of the probe is measured and compared with the 
input state. From the comparison, if we know the physical 
mechanism that governs the combined probe-target dynamics 
(e.g. the interaction Hamiltonian), we can deduce the value of 
the parameter. In general, the measurement process is affected 
by statistical errors, whose origin can be extrinsic (e.g. envi¬ 
ronmental noise) or intrinsic (e.g. Heisenberg uncertainty rela¬ 
tions, input and output states being in general non-orthogonal 
and hence not distinguishable with certainty). 

To improve the precision of the estimation, several strate¬ 
gies can be adopted. First, we can optimize the input state of 
the probe so that the probe-target interaction is able to imprint 
the highest possible amount of information about the target 
parameter into the probe, i.e. the input and output states be¬ 
come most distinguishable. In particular, there might be states 
of the probe that are left unchanged by the interaction with the 
measured system and are useless in this sense, so we usually 
want to avoid them. Second, we can repeat the measurement 
several times to enlarge our statistical ensemble of data and 
extract a sharper expectation value. This can be realized by 
preparing many copies of the probe and making them interact 
independently with the system (parallel scheme), or by mak¬ 
ing the same probe interact repeatedly with the system before 
extracting the information (sequential scheme). Third, we can 
exploit the presence of genuine quantum resources, such as 
quantum coherence, or quantum correlations either between 
the many copies of the probe or between the probe and some 


ancillary system that is kept as a reference, to gain advantage 
over purely classical strategies. In particular, it is well known 
that the presence of entanglement allows one to estimate a 
parameter encoded in a unitary dynamics (e.g. a phase shift) 
with an error that scales as IjN with respect to the number 
N of collected measurements, while classical strategies can at 
most achieve a scaling of 1/ vJn [3, 5]. 

In some specific cases of practical relevance, we may not 
have a complete trusted knowledge of the probe-target inter¬ 
action mechanism and therefore we may find it harder to op¬ 
timize the input state of the probe in order to maximize the 
efficiency of the estimation. For example, we could imag¬ 
ine a situation in which we become aware of unwanted noise 
sources just before we retrieve the output state, meaning that 
the actual transformation is different from what we expected 
when we prepared the probe, which is then likely to be sub¬ 
optimized. As another example, we could be asked to pre¬ 
pare a passe-partout probe state that must be good whenever 
the interaction with the measured system is described by a 
Hamiltonian picked at random from a given ensemble, so that 
we have no interest in optimizing the probe for a particular 
element of the ensemble. It turns out that in such and sim¬ 
ilar situations, that we may describe as instances of “black¬ 
box” quantum metrology, the presence of correlations gives 
another fundamental advantage [6-10]. While with a single 
probe system we always run the risk of preparing the probe in 
a state which is left unmodified by some unlucky interaction 
mechanism with the target system, by exploiting correlations 
between the probe and an ancillary system kept as a reference 
we can instead guarantee a minimum detection efficiency. 

It is then interesting to ask the following question; Given a 
certain minimum efficiency that we want to achieve in a black¬ 
box quantum metrology task, what resource should we look 
for in our probe state? The answer has been found in several 
recent works [6-10] and in short is: discord-type correlations. 
These are general quantum correlations that encompass en¬ 
tanglement but also describe the nonclassical nature of most 
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separable states. They have been introduced for the first time 
in 2001 under the name of quantum discord [11, 12] and have 
been the subject of extensive studies in the last decade [13]. 
In particular, it has been recently shown that quantum cor¬ 
relations in a bipartite probe can be exploited to guarantee a 
minimum precision in the estimation of a local phase [7, 8] 
or a minimum probability of detecting a remote object in a 
quantum illumination [9] or quantum reading [10] scenario. 
Let us stress the following fact. While, as one could expect, 
pure maximally entangled states of the probe-ancilla bipartite 
system are still the best option for the considered tasks, en¬ 
tanglement is not a necessary resource in the black-box sce¬ 
nario. On the contrary, discord-type correlations embody the 
fundamental feature that provides, guarantees and quantifies a 
quantum over classical advantage in a vast class of metrology 
tasks (see also [14]). Therefore, one can also consider using 
“cheaper” separable but quantumly-correlated states [15, 16] 
if the required minimum precision is not too stringent, and in 
general if the production of pure entangled states is hindered 
by technological limitations. 

In this paper we extend the above analysis a significant step 
further. As just discussed, the amount of discord-type quan¬ 
tum correlations in the input state of the probe is all the in¬ 
formation that we need in order to know what the worst-case 
performance will be and hence guarantees a minimum esti¬ 
mation efficiency. However two states with the same amount 
of discord-type correlations are not fully equivalent resources 
from a general metrological point of view. Indeed, although 
they are characterized by the same minimum estimation ef¬ 
ficiency, one of the two states could be better on average 
and thus preferable over the other, as long as the informa¬ 
tion about the system-target interaction remains partially un¬ 
known. For all practical purposes, truly versatile probes for 
quantum metrology should then be able to offer acceptable 
performances on average when employed for a broad range 
of tasks. Therefore, other than investigating the resources in¬ 
volved in determining a worst-case performance as done ear¬ 
lier, one should address a different key question; Given a cer¬ 
tain average efficiency that we want to achieve in a black-box 
quantum metrology task, what resource should we look for 
in our probe state? Here, we discuss this aspect in full de¬ 
tail and we provide a comprehensive classification and char¬ 
acterization of bipartite quantum probe states in terms of their 
average metrological performance. Together with previous re¬ 
sults [6-10], our analysis can have a direct impact on the con¬ 
crete search for optimal and versatile probe states useful for a 
plethora of metrological applications in realistic conditions. 

To deliver a quantitative analysis, we focus here on the skew 
information I(p, H) - -Tr |^[ ^Jp, //]^j /2, which expresses the 
amount of information stored in a state p that cannot be ac¬ 
cessed by measuring the observable H, due to the noncommu¬ 
tativity between state and observable [17, 18]. The skew infor¬ 
mation is one possible extension of the classical Fisher infor¬ 
mation to the quantum domain, being part of a larger family 
of Riemannian contractive metrics on the quantum state space 
[19, 20]: therefore, it directly quantifies the susceptibility of a 
probe state p to an infinitesimal change in a target parameter 
encoded in the observable H. If the observable acts locally 


on one subsystem of a bipartite state, the skew information is 
bounded from below by the amount of discord-type correla¬ 
tions in the state and its minimum value can be used in fact 
as a measure of discord-type correlations, defined in [6] as the 
Local Quantum Uncertainty (LQU). This quantity is closely 
related to other measures, such as the Interferometric Power 
(IP) [7] and the Discriminating Strength (DS) [9], that have a 
direct interpretation in terms of metrological tasks in worst- 
case scenarios. For example, the LQU coincides with the DS 
for qubit systems and gives a lower bound to the IP in gen¬ 
eral. Therefore the LQU can be interpreted as a minimum 
susceptibility of a bipartite state to local transformations on 
one subsystem, thus being relevant from a quantum estima¬ 
tion perspective. Moreover it is based on a simple functional, 
the skew information, that is typically easy to compute and 
serves as a good starting point for our investigation. 

For arbitrary states of a generic bipartite system, we com¬ 
pute here the average of the skew information over specific 
classes of local observables acting on one subsystem. The 
resulting quantity, referred to as Average Skew Information 
(AvSk), quantifies therefore the average susceptibility of a bi¬ 
partite state to local transformations. Remarkably, such an av¬ 
erage susceptibility can be expressed through a simple analyt¬ 
ical expression, that clearly shows what is the role played by 
the properties of the observables and by the properties of the 
state in determining the average performance. Thanks to this, 
we provide an extensive characterization of the AvSk and of 
its features. In the specific case of a two-qubit system, where 
the LQU is also computable in closed form [6], we then carry 
out a parallel study of our new quantity and of the LQU that 
allows us to identify which states of the probe are better given 
different constraints. It turns out that the resources needed 
in the probe state to optimize the average metrological perfor¬ 
mance are quite distinct from those (discord-type correlations) 
needed instead to guarantee a minimum performance. We also 
find that our AvSk is equivalent, up to a numerical prefactor, 
to another quantity recently introduced by Luo and cowork¬ 
ers [2 1 ] which is similarly based on the skew information but 
considers a different kind of averaging. This connection al¬ 
lows us to easily prove that the AvSk can be adapted to define 
a measure of correlations but not specifically of quantum (like 
the LQU) or classical correlations. Furthermore, our analy¬ 
sis complements that of Luo et al. by finding a nice closed 
analytic expression and a clear operational meaning for their 
measure. Finally, we also compute the variance of the skew 
information to investigate what additional knowledge can be 
gained from higher moments of the statistics. 

The main content of the paper is structured as follows. In 
Sec. II we compute the average of the skew information over 
an ensemble of local observables with fixed non-degenerate 
spectrum. In Sec. Ill we enumerate and prove the basic prop¬ 
erties of the average skew information. In Sec. IV we discuss 
how the average skew information depends on the choice of 
the spectrum of the local observable. In Sec. V we compute 
the average skew information for specific classes of states and 
we derive some general bounds. In Sec. VI we make a detailed 
analysis of the two-qubit case, comparing the average skew 
information with the LQU (i.e. the minimum skew informa- 
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tion). In Sec. VII we also compute the variance of the skew 
information and we discuss what this refined statistics can tell 
us about the presence of quantum correlations. In Sec. VIII 
we discuss the connection between our quantity and the one 
recently introduced by Luo et al. [21], and we provide addi¬ 
tional comments on the role of correlations. Finally, in Sec. IX 
we provide an explicit interpretation of the main results of 
this paper from a metrological point of view. We present our 
concluding remarks in Sec. X. Some technical derivations are 
deferred to Appendices. 

II. AVERAGE OE THE SKEW INEORMATION OVER 
LOCAL OBSERVABLES WITH FIXED NON-DEGENERATE 
SPECTRUM 

If p is a density operator on a Hilbert space 'Kz and H is 
an Hermitian operator on 'Hx, the skew information of p with 
respect to H is defined as [17, 18] 

( 1 ) 

and expresses the amount of information stored in a state p 
that cannot be accessed by measuring the observable H, due to 
the noncommutativity between state and observable. Note that 
in general it is always possible to find an observable Hp which 
is diagonal in the eigenbasis of p and therefore can grant com¬ 
plete knowledge of the state, i.e. I{p, Hp) - 0. However, this 
is no longer true if we make the additional assumption that 
observables act only on a part of the global system. 

It has been recently shown [ 6 ] that when p - p^s is a den¬ 
sity operator of a bipartite system described by the Hilbert 
space HIab - ® and = /Li ® is a local Hermitian 

operator acting only on HI a, the skew information is bounded 
from below by the presence of general nonclassical correla¬ 
tions of the discord type [11-13] in the state p. Quantum dis¬ 
cord, as proposed in the original formulation [ 11 , 12 ], mea¬ 
sures the part of the information stored in the correlations of 
a bipartite system AB that cannot be retrieved by measuring 
locally one of the subsystems (say A). This locally unaccessi- 
ble correlations arise because a local measurement can perturb 
the state of the system by projecting it onto a particular local 
basis for A, losing some information in the process, and the 
existence of an unperturbing measurement is not guaranteed. 
In the same spirit, taking the minimum of the skew informa¬ 
tion over some ensemble of local observables of a bipartite 
system gives the minimum incompatibility between the state 
p and the ensemble of observables, i.e. the amount of infor¬ 
mation that always remains hidden under a certain family of 
local measurements. In particular, if one considers the set of 
all local observables with a fixed non-degenerate spectrum, 
one obtains the LQU introduced in [ 6 ] 

min I{p,H(^a)) , ( 2 ) 

|H(Aa)) 

which is in fact a good quantifier of discord-type correlations. 
In Eq. (2) the minimum is taken over a set of local observ¬ 
ables with fixed non-degenerate spectrum A .4 = 


where (|/)a) is an orthonormal basis of A and the /l,’s are 
all different. This is necessary to ensure that the identity Ia 
is excluded from the minimization set and the trivial case 
I(P,1a) = 0 is avoided (this must hold also if considering 
any subspace of HIa)- That is, only observables of the form 
H(Aa) = UaAaUa^ are considered, where Ua is any local 
unitary transformation on subsystem A. As shown in [6], the 
LQU satisfies all the properties required to a well-behaved 
measure of discord-type quantum correlations [22, 23]. In 
particular it is zero if and only if the original quantum dis¬ 
cord is zero and hence captures the same type of correlations. 
Moreover, the LQU is strongly connected to other measures 
of quantum correlations, such as the IP [7] and the DS [9], 
that have a clear interpretation in a metrological context. For 
example, the LQU coincides with the DS if the bipartite sys¬ 
tem is made of two qubits, and in this case it measures the 
minimum efficiency of a given bipartite state as a probe for a 
quantum illumination task [24] where one must decide if any 
transformation in a given set of isospectral local unitary oper¬ 
ations has been performed or not on the probe. 

Here, instead of taking the minimum as in Eq. (2), we com¬ 
pute the average of the skew information over the set of Her¬ 
mitian operators UaAaUa' spanned by the unitary group on 
Hi A- In light of the above discussion, this quantity, which will 
be named simply Average Skew Information (AvSk), can be in¬ 
terpreted as the average susceptibility of a bipartite probe to 
local transformations and local parameters. The AvSk can be 
written as an integral with respect to the Haar measure of the 
unitary group dpH(UA) 

I^^ip)^ J dpH{UA)I(p,UAAAUA^) 

= - ^ J dpH{UA)^r [[ Vp, UaAaUa^]"] . (3) 

In choosing our notation, we made explicit the fact that the 
AvSk depends only on the state and on the specific choice of 
the spectrum. To compute the integral in Eq. (3) we start by 
rewriting Eq. (1) for the case of a bipartite state p - pab and a 
local observable 7/ = Ha ® Ib as 

/(p, Ha) = Tr[( yffHA){HA yff) - ( #Ha)( >/^//a)] 

= Tr[( yf^HA ® Ha- y/pA'B' 

- ^PabHa ® ylpA’B’HA’)S ab\a’B’], ( 4 ) 

where following the procedure of Ref. [25] we introduced 
a copy HIa'B' - 'Ha' ^^Hb' of the original Hilbert space 
HIab -'Ha®'Hb and the swap operator 5 ab|A'B' acting on 
Hab ®'Ka'B' [26]. Using Eq. (4) and the properties of the 
swap operator (see Appendix A) we can now rewrite Eq. (3) 
as 

J^^(p) = Tr[(pAB ® Ia'B' - nW ® VPa'B' ) 

x7-®(Aa®Aa-)5ab|a'B']. (5) 

where 7 “*^^(Aa ® Aa-) is the so-called twirling channel [27- 
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29] applied to the operator Aa ® (see Appendix B) 

= J dfiH(UA)(UA ® I/aOCAa ® Aa>)(uI ® t/],) 

_ AATr[A^]^-Tr[A^] ^ 

Na(NI - 1) 

AATr[A 2 ] - Tr[A ^]2 

" mK - 1 ) ® 

In writing Eq. ( 6 ) we introduced the dimension Na of the 
Hilbert space 'Ha- Plugging the last two lines of Eq. ( 6 ) into 
Eq. (5), using again the properties of the swap operator, and 
evaluating the trace, we finally get a remarkably compact for¬ 
mula for the AvSk of an arbitrary bipartite state p. 


jA.(^)^^^Tr[A 2 ]-Tr[AJ 


NAiNi - 1 ) 


2 r 

- A,4-Tr4(Tr4[^/^])2] 


(7) 


We stress that the analytic expression Eq. (7) holds for any 
dimension of the Hilbert spaces 'Ha and 'Hb- 


III. PROPERTIES OE THE AVERAGE SKEW 
INFORMATION 


We discuss now some properties of the AvSk 

Property la - Eor any fixed spectrum, the AvSk is 
non-negative. This is trivially true as the skew information is 
non-negative and this is not changed by taking the average. 

Property lb - Eor any fixed non-degenerate spectrum, the 
AvSk is zero if and only if the state is of the form 

I/t 

Pab - — ®Pb- ( 8 ) 

A a 

The proof of this is rather long and is postponed to Sec. V. 

Property 2 - The AvSk is invariant under local uni¬ 
tary operations Wa, Vb- Indeed, consider the transforma¬ 
tion p —» {Wa ® ® Eg) which also maps yfp into 

{Wa ® Vb) Then, by exploiting the cyclic prop¬ 

erty of the trace in Eq. (7), it is easy to see that 

{{Wa ® Vb)P (W] ® V^)) = (p). (9) 


Property 3 - The AvSk is non-increasing over all com¬ 
pletely positive trace-preserving (CPTP) maps acting locally 
on B. To show this, let us first decompose an arbitrary local 
CPTP map < 1 >b(p. 4 b) as a unitary interaction with an external 
environment followed by a partial trace over the degrees of 
freedom of the environment [30] 

®b(Pab) = Ti'b {Ube(Pab ® PeW^be] ’ 


where we can further assume that the unitary operation in¬ 
volves only the subsystem B and the environment, without af¬ 
fecting A. Mimicking the demonstration of Property 2, we can 
show that the skew information satisfies the following prop¬ 
erty 


I(Ube(Pab ®Pe)U\^,Ha) — I{pAB ®Pe,Ha), (11) 

and it is also easy to see that 

KPAB ®Pe, Ha) = I{pAB, Ha). (12) 

Einally, it was proven in [31, 32] that 

I{^b{Pab),Ha) = ^ {T'J^e[Ube{Pab ® pE)Ug^^, Ha) 

^ I{Ube{Pab ® Pe)U\^, Ha) 
^I{Pab,Ha). (13) 

Since Property 3 is true for the skew information itself, it 
remains true also when taking the average. 


Property 4 - Eor pure states, the AvSk is an entanglement 
monotone. Indeed, given a pure state \i//}ab, we have that 
VI'A)ab<iAI = I</')ab<<AI- Plugging this into Eq. (7) leaves us 
with 


I'^nmAB) = 


AATr[A^] - Tt[Aa? 
Na{NI - 1) 


Na - Trafpl] 


(14) 


where pB - Tila[\4>)ab{4'W- For pure states, a convenient mea¬ 
sure of entanglement is provided by the generalized con¬ 
currence C{\\}/)ab) [33], which depends only on the pu¬ 
rity of the marginal density operators TrB[p|] = TrA[p^] as 

C{\4>)ab) - .^2 - 2TrB[p|]. We can then rewrite 


I^H^ab) = 


AATr[A^]-Tr[AJ^ 
Na{NI - 1 ) 


C{W)ab)^ 

2 


(15) 

which clearly makes the AvSk an entanglement monotone. 
Note that, however, it cannot be considered strictly speaking 
as a measure of entanglement since it does not vanish on all 
separable (product) states; still, one can obtain a fully fledged 
entanglement measure on pure states by rescaling the AvSk 
subtracting the dimension-dependent constant in Eq. (15). 


Property 5 - The AvSk is convex with respect to the state. 
The result follows simply from the convexity of the skew in¬ 
formation [17], that is preserved by taking the average, 

I^Hpp, + (1 - p)pi) < P I^Hpi) + (1 - p) I^Hpi) , (16) 

for any two states p\,p 2 and a probability 0 < p < 1. This is 
noteworthy, since convexity is lost instead when taking the 
minimum rather than the average [ 6 ]. 


Prom Property 1, we immediately see that, at variance with 
the LQU, the AvSk is not a proper measure of discord-type 
quantum correlations [ 22 ], because in general it is different 
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from zero when evaluated on the set of classical-quantum or 
classical-classical states [13]. In fact, the AvSk can be non¬ 
zero even for completely uncotTelated states, e.g. for any state 
of the form \ip)A ® \^)b 7 meaning that it is neither a measure 
of classical nor total correlations. Nevertheless, we will see in 
Sec. VIIIB how to construct a proper measure of total corre¬ 
lations based on a modification of the AvSk, recovering and 
complementing the analysis done in Ref. [21]. We remark 
however that our focus here is not to define yet another ab¬ 
stract measure of correlations. Instead, we are going to use 
the AvSk operationally as a guidance to identify optimal probe 
states for (black-box) quantum metrology, adopting their av¬ 
erage performance as our figure of merit. 

In the following Sections we are going to compute the AvSk 
for some specific classes of states and derive some general 
bounds, that can be straightforwardly established thanks to 
convexity. Finally, if we study simultaneously the AvSk and 
the LQU, we can point out what states are better used in quan¬ 
tum sensing and metrology tasks such as state discrimination 
and parameter estimation, depending on the rules of the game. 


IV. DEPENDENCE OE THE AVERAGE S KE W 
INEORMATION ON THE SPECTRUM 


The expression (7) for the AvSk that we found at the end 
of Sec. II, explicitly factors the dependence on the spectrum 
of the observable and the dependence on the state. In this 
Section, we investigate how different choices of the spectrum 
relate to one another. 


A. Invariance under translation of the spectrum 

First we show that if two spectra A^ and A^ are connected 
by a rigid shift, the two induced AvSks are equal. The rigid 
shift condition is expressed as A^ = Aa -i- where 77 is any 
real number. We then have 

Tr[A;,] = TiIAa] + tiNa, (17) 

Tr[A'/] = Tr[Al] + 277Tr[AA] + ii^Na. (18) 

Plugging the above expressions into Eq. (7) and considering 
only the part containing the spectrum, we easily see that 

jVATr[A;^]-Tr[A;]^ ^ NaTi^Aa^] - Tr[AA]^ 

Na(NI-1) Na{NI-1) ■ ^ ^ 

This implies that V77. Therefore, this al¬ 

lows us to simplify Eq. (7) by considering only spectra with 
trace equal to zero. 


TiIAa] = 0 


J^^(p) - 


Tr[A^] 

Nj-l 


VA-Tr4(TrA[>/^])2] . 

( 20 ) 


B. Scaling under scalar multiplication of the spectrum 


Next, we consider what happens if we take a spectrum Aa 
and transform it to tjAa by scalar multiplication. Thanks to 
Eq. (20), it is immediate to see that - if for 

any value of 77 . 


C. Optimal spectrum 


We can now ask which spectrum yields the highest prefac¬ 
tor to the AvSk. Erom the previous results, it is obvious that 
multiplication of a spectrum by a big real number can make 
the prefactor as big as desired. However, we want here to 
highlight the role played by the distribution of the eigenvalues, 
rather than their magnitude. We can make a fair comparison 
by exploiting the translation invariance and the scaling intro¬ 
duced above, and considering only positive spectra with unit 
trace (i.e. we map each spectrum to a density matrix). We then 
see from Eq. (7) that all the information about the spectrum is 

in the prefactor , which for a fixed dimension Na de- 

pends only on the spectrum purity. Therefore it is immediate 
to see that the best spectra are those that have Na - 1 degen¬ 
erate eigenvalues, i.e. those spectra that can be mapped into 
pure-state density matrices by means of rigid shifts and scalar 
multiplications. Eor example one such spectrum, taken trace¬ 
less to satisfy the condition discussed in Sec. IV A, is given 
by Aa = ((Na - 1)/Na, -1/Na, ..., -1/Na]. 

This means that if we want to encode some information on a 
state but we cannot choose the encoding basis, an almost fully 
degenerate spectrum allows to encode, on average, the max¬ 
imum amount of information. We stress that this situation is 
almost opposite to what happens for the LQU [ 6 ] and for simi¬ 
lar measures of quantum correlations such as the IP [7] and the 
DS [9] that consider the worst-case performance, where it is 
instead believed that the optimal spectrum is harmonic [9, 34], 
i.e. fully non-degenerate and with equally spaced eigenvalues. 
Eurthermore, we see that the AvSk is non-trivial as soon as 
the spectrum has some different eigenvalues, i.e. as soon as 
Aa 4^ Ia. We don’t need to impose here the stricter condition 
of full non-degeneracy required, for example, by the LQU. 


V. DEPENDENCE OF THE AVERAGE SKEW 
INFORMATION ON THE STATE 


In this Section, we study the AvSk for specific classes of 
states or, conversely, we look for the states that yield the max¬ 
imum and the minimum AvSk given specific constraints. All 
the results provided here hold for any dimension of TYa ® TYb. 
Without loss of generality, we consider traceless spectra (see 
Sec. IV A). 
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A. Average skew information for pure states 

We start by considering pure bipartite states. As we have 
seen in Sec. Ill, the AvSk takes a simple form on the set of 
pure states \4i)ab. 




Tr[A2] 


Nl-l 


Na - Trglpl] 


( 21 ) 


where pg - Txa[\4')ab{4'W is the reduced state of subsystem 
B and Tr^lpg] is its purity, which can take values between 
1/ mm{NA,Ng} and 1. Therefore we can find the following 
bounds for the AvSk of pure states: 


Tr[A2]. , . Tr[A2] 

- l] < < 


N- 


Ni-l 


Na- 


1 


min{AA,AB) 

( 22 ) 

where the upper bound is saturated by pure maximally entan¬ 
gled states and the lower bound is saturated by pure product 
states. 


B. Average skew information for separable states 

Another interesting class of states is given by separable 
states. Here we have no entanglement and we can investi¬ 
gate if the presence of discord-type quantum correlations has 
a specific impact on the AvSk, as it has for the LQU [6]. We 
start by considering a general separable state 

Psep = ^P/Pa (23) 


where p® and p® are arbitrary density matrices of A and B, 
Pi > 0 and 2, Pi - 1. From the convexity of the AvSk (see 
Sec. Ill), we have 

0 < J^-'fpsep) < y Pi -^"^^(pa ® Pb ) ^ max I^^ipA ®pg) , 

' ' ■ A -' ' ' {pa^pb} 

(24) 


product states where the local density matrix on A is pure, i.e. 

p = I<A>a<iAI ®pb- 

A few remarks are in order here. First of all, we notice that 
all separable states yield a lower AvSk than any pure entan¬ 
gled state. We can then use the AvSk as a witness of entangle¬ 
ment and say that 

TrlA^ 1 

I^^(p) > — —^[Aa - l] p is entangled. (27) 

Furthermore, since the maximum AvSk among separable 
states is reached by a completely uncorrelated state, we can 
claim that the presence of quantum correlations other than en¬ 
tanglement has no specific effect on the average susceptibility 
of a bipartite state to local transformations. Of great impor¬ 
tance is instead the local purity of the probing subsystem A: 
as soon as pA is not maximally mixed, an average metrologi¬ 
cal performance is guaranteed even in absence of a correlated 
reference subsystem B. 

We recall, however, that discord-type correlations as mea¬ 
sured by the LQU determine instead the minimum suscepti¬ 
bility of a bipartite state to local transformations. A compar¬ 
ative analysis of the AvSk and of the LQU can serve then to 
identify states that simultaneously yield satisfactory levels of 
complementary figure of merits and emerge as suitable probes 
for sensing applications. We will come back to this point in 
Sec. VI, where we investigate the specific case of two qubits. 


Average skew information for classical-quantum states 


We compute here the AvSk for a specific class of separable 
states, i.e classically correlated states that have zero LQU (or 
equivalently zero quantum discord). Since we are considering 
local measurements on subsystem A, the set of classically cor¬ 
related states is given by the so called classical-quantum (CQ) 
states [11, 13] 


where in the last term we take the maximum over all product 
states p - Pa® pg. By direct substitution in Eq. (20), we have 

A Tr[A2]|- ,, 

0 < I^Hpa ®pg) = Aa - Tr[ 

A;^ - 1 J 

Tr[AA]r i 

^ - 4 ( 25 ) 

and finally 

0 < (psep) < [Aa - l]. (26) 

The lower bound is saturated, for example, by product states 
of the form p = Ia/Aa ® pg (as announced in Sec. Ill, and 
as we are going to show, these are the only states with zero 
AvSk) while the upper bound is saturated, for example, by 


Na 

PCQ = Pi IOa <*1 ® (28) 

1=1 

where (p,) is a set of probabilities, (IOa) is an orthonormal ba¬ 
sis of A and jp®) are general density matrices for subsystem 
B. Note that for any such state the existence of a commut¬ 
ing local observable Ha that nullifies the skew information is 
guaranteed (i.e. when Ha is diagonal in the basis (IOa))- The 
CQ states include the so called classical-classical (CC) states 

Na Nb 

Pcc - ^29) 

1=1 j=l 

where now also {ly)^) is an orthonormal basis of B. Starting 
from an arbitrary CQ state, we plug Eq. (28) into Eq. (20) and 
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get 


Tr« 


(TtA [ -v^]) 


= Tr 


Na Na 

Yj P‘P^b + 2 Y^Pi ^P'b iPB 

j>i=i 


i=l 


Na 


^1+2 Y yfpip^i Tr 

i>i=l 


kwpi" 


(30) 


A lower bound to Eq. (30) is given by 
\2t 


Tr« 


(Tr4^]) 


= 1 . 


The bound is saturated, for 


example, when only one of the /7,’s is non zero, i.e. for 
product states (i/rl ® ps- Another possibility is that the 

.(0 


p^’s are all orthogonal to each other. For example, the set 

(p®) could be a set of pure orthogonal states {\([>i)g (0,1) on B 
(thus giving a CC state). The corresponding upper bound to 
the AvSk of CQ states becomes 




Tr[A^] 
Nl- 


-1) 


(31) 


We can also hnd an upper bound to Eq. (30) if we use the 

\21 

> 0, namely 




inequality Tr 

1 + 2 ^ VW Tr 

;>/'=! 

Na 

<1+ £ V]^(Tr[p®] + Trp> 

j>i=i 

Na Na 

= 1+2 v^= Y 


Na 


u=i 


j>i=i 


Pi + Pj 




^Na. 


(32) 


where {|/)b) is an orthonormal basis of B and jp®) are general 
density matrices for subsystem A. We plug this into Eq. (20) 
and get 


1 <TrB 


(Tr4#^]f] = f;p,(Tr 

i ' 



(35) 


The lower bound is saturated if and only if all the (p®) in 
Eq. (34) are pure states, i.e. for all density matrices that can 
be written in the form 


Na 

PpQC = Y P' 

1=1 

where is a set of generic pure states of A (in particu¬ 

lar we don’t require them to be orthogonal, at difference with 
the set (IOb)). We will use the name (pure quantum)-classical 
(pQC) for states of the form Eq. (36). We stress that PpQc is 
not itself pure in general (that’s why we put the word “pure” 
between parenthesis in the full name and we write a small “p” 
in the abbreviation). The upper bound is saturated if and only 
if all the (p®) are proportional to the identity, i.e. again for 
states of the form Ia/Na ® Pb- Correspondingly, for the AvSk 
we get 


0< J^^(PQC)< ^£^[a^a-1]. (37) 

As anticipated, we see that the AvSk on the subset of QC states 
achieves the same bounds as the AvSk on the set of CC states. 
This means that general quantum correlations have no clear 
effect on the average susceptibility of the state. Instead, we 
see again that the purity of the local state of A has a great 
importance. 


C. Maximum and minimum of the average skew information 
for general states 


The bound is saturated if and only if p, = 1 /Na for each i and 
all the p^^’s are equal. In other words, the CQ state must be 
of the form Ia/Na ® Pb to have zero AvSk. In conclusion, the 
bounds to the AvSk of CQ states become 

0< J^^(Pcq)< ^^[iVA-l]. (33) 


Average skew information for quantum-classical states 

We can also compute the AvSk on the set of quantum- 
classical (QC) states that, opposite in spirit to the CQ states, 
can have a hnite amount of discord-type quantum correlations 
(as measured e.g. by the LQU with respect to subsystem A). 
An arbitrary QC state can be written as 

Na 

PQC = Y P‘P^A ® ’ (34) 

1=1 


Another interesting question is what states, including po¬ 
tentially entangled states, have the absolute highest and low¬ 
est AvSk. Consider a general bipartite state p and its diagonal 
expansion onto some basis {|i/'i)ab)- 


p^YPi l‘/'i)AB (<Ail, —» I'Ai)ab (^i\ ■ 

(38) 

From the convexity of the AvSk we have 


Pi\>/'i}ABm) < Y P‘ I'^^mABm 


N^a- 


l[«'> 


1 


mm{NA,NB 


], (39) 


where the absolute maximum of the AvSk for pure states, and 
hence for all states, is reached only by the maximally entan¬ 
gled states. 

We look now for the minimum. From the very dehnitions 
of the LQU and of the AvSk, we have the simple relation 
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Therefore, the states with minimum AvSk 
can only be found within the set of states with minimum (zero) 
LQU, i.e., the CQ states. As we have already seen, among all 
the CQ states only the states of the form IaINa ®Pb have zero 
AvSk. This gives a proof of property lb, that we formulated 
in Sec. III. 


D. Minimum of the average skew information for fixed LQU 

Another interesting question is what states have AvSk equal 
to their LQU. We recall that the LQU can be expressed as 


We compute the AvSk and the LQU (using the formula in [6]) 
for 10^ randomly generated two-qubit states. In Fig. 1 we plot 
the AvSk of each state vs. the corresponding LQU. 

The results of Sec. V are clearly illustrated by the plot. 
Namely, we observe the following: 

• Since the LQU is obtained through a minimization over 
all possible unitaries and the AvSk is obtained through 
an average, we must have that I°'^(p) > 'l(°'’(p). This 
lower bound, shown by a blue solid line in Fig. 1, is 
saturated, for example, by isotropic and Werner states 
(see Appendix C). 


(p) = 1 - Tr [ ^f^HA ^f^HA ], (40) 

where Ha is some Hamiltonian minimizing the skew informa¬ 
tion. The AvSk over all Hamiltonians with the same spectrum 
can be equal to the LQU if and only if 

Tr [ ^Ha yf^HA] = Tr [ ^{uIHaUa) ^[^{uIHaUa)] 

= Tr [(Ua ^/^uI)Ha(Ua ^uI)Ha] , VU^. 

Note that we can always add an arbitrary local unitary trans¬ 
formation Vb to the density matrix without affecting the above 
equality. In other words, we can ask that 

Tr [ ^/pHa ^/pHA] = Tr [(Ua ® Vb yfpU[® v]^)Ha 

X (Ua ® Ub Vp vI)Ha\ , Wa- (41) 

A sufficient condition for Eq. (41) is expressed as 

VUa, 3Ub such that s/p - (Ua ® Vb sfp U\® V^). (42) 

Therefore, we must look for states that are invariant un¬ 
der any unitary operation Ua if we allow the application 
of an arbitrary local correction Vg- Some examples, when 
Na - Nb, are given by the Werner states [35], that sat¬ 
isfy pw - (U ^ U)pw(U ® U)'' Vt/, and by the isotropic 
states [27], that satisfy pi - (U ® U*)pi(U ® t/*)^ Vt/. 


• The separable states, including the CQ states (for which 

the LQU vanishes) and the QC states, satisfy the bound 
{j-^(p,ep),J°'=(pcQ),J‘"'(pQc)} < 2/3. CQ states 
are shown by a green dot-dot-dashed line in Fig. 1. 
pQC states have all = 2/3 and are shown by 

the horizontal dashed black line. CC states all have 
'^/‘"^(pcc) = 0 and I'^^(pcc) = 2/3. 

• From Ref. [9], we know that the separable states must 
have limited LQU, 'Z/°’'(psep) < 1/2. Therefore, they 
must lie left of the vertical dashed black line in Fig. 1 . 
Combined with the previous observation, this allows us 
to identify a region where only entangled states exist 
and a region where separable and entangled states co¬ 
exist (see also appendix C). 


• The pure states satisfy the bound 2/3 < I°''^(\4i)ab) ^ 1- 
The lower bound is saturated by separable pure states 
and the upper bound is saturated by maximally entan¬ 
gled states (Bell states). The Bell states also achieve the 
highest AvSk among all states. Moreover, for any given 
value of the LQU, the highest possible AvSk is achieved 
by a pure state. For pure states of two qubits, we have 


I‘^^(\/j)ab) - 3 


1 -H 


^‘^-'(\/f)AB) 


C(W)ABy 


, (44) 


VI. AVERAGE SKEW INFORMATION FOR TWO QUBITS 

We focus now on the exemplary case of two qubits, for 
which the analysis becomes particularly simple and insight¬ 
ful. Indeed, in this case we can also explicitly compute the 
LQU [6] and we can classify all the states according to their 
minimum and average susceptibility to local transformations, 
looking at the results of Sec. V in more detail. Furthermore, 
the LQU of two qubits coincides with their DS [9], and the 
AvSk can be then rigorously interpreted as the average dis¬ 
crimination efficiency of the state in a quantum illumination 
task [24]. The analysis takes then an explicit metrological 
connotation. 

From the results of Sec. IV we can fix A^ = cr^ without 
loss of generality, where cr^ is the third Pauli matrix, and the 
expression of the AvSk for any two-qubit state becomes then 


where C(|i/')ab) is the concurrence. Pure states are indi¬ 
cated by a red dot-dashed line in Fig. 1 . 

• States of the form (Ia/2) ® pB are the only states having 
zero AvSk. 

The simultaneous analysis of the AvSk and of the LQU pro¬ 
vides a useful guide when we need to decide which states 
of the two-qubit probe are more suitable to perform a given 
metrological task (e.g., in the present case, state discrimina¬ 
tion). We immediately see that maximally entangled states, as 
can be expected, are the best choice when we focus on both 
the worst-case performance and the average performance as 
figures of merit. However, if we have limited resources and 
do not have access to entangled states, we can still achieve 
good results using separable states. For example, the state 


r^^ip) = 


TrB[(Tr^[#])2] 


(43) 


P = ^| 0 >A< 0 | ® | 0 )b< 0 | + ^|+)a<+I ® IDbQI ( 45 ) 
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Figure 1. (Color online) AvSk and LQU for 10^ randomly generated two-qubit states (gray dots). Special classes of states are highlighted by 
different lines and detailed in the sketch on the right. See also the main text for a complete description of the various regions and boundaries. 


yields a LQU equal to 1 /2 and an AvSk equal to 2/3 (recall 
that the maximum is 1 for both quantities). This state, among 
all separable states, has the highest amount of discord-type 
correlations. This confirms that quantum correlations beyond 
entanglement are indeed useful for metrological applications, 
although they play a relevant role only in determining the 
worst-case performance but have little effect on the average 
performance (the value 2/3 for the AvSk can be reached even 
with product states). Another observation that we can make 
is the following. If one needs to guarantee a minimum effi¬ 
ciency of the probe, i.e. fix the LQU as a primary figure of 
merit, there is still some freedom in the choice of the initial 
state, with pure states being on average better than any other 
possibility. Our analysis of the AvSk can be very useful in this 
sense. 


VII. VARIANCE OF THE SKEW INFORMATION 

In this Section, we complement the above analysis by com¬ 
puting the variance of the skew information, which tells us 
how much the efficiency of a given probe state is fluctuating 
around the average value for different choices of the encoding 
unitary. The variance is defined as 

AJ^-(p) = J diUH(UA)f(p, UaAaUa^') - .(46) 

The first term on the right-hand side of Eq. (46) can be com¬ 
puted following the prescriptions of Appendices D and E, and 
the expression for two-qubit states is given by Eq. (El3). We 
can further impose that A^ = cr^ without loss of generality. 
Since analytic insight is out of reach for such a cumbersome 
expression, we resort again to computing the variance numer¬ 
ically for the 10^ random two-qubit states generated before. 

The results are presented in Eig. 2. On the left panel, we 
show a density plot of the square root of the variance, given 


the corresponding AvSk and LQU. It is immediate to see that 
the variance is zero on the isotropic and Werner states, for 
which the average is equal to the minimum, and is bigger 
when the difference between the average (AvSk) and the min¬ 
imum (LQU) is bigger, as should be expected. Moreover, we 
find that there are precise quantitative relations that describe 
this behavior. To show this, on the right panel of Eig. 2 we plot 
the square root of the variance versus the corresponding value 
of the difference (AvSk - LQU). All the points lie within a 
well-defined region and we can find the states lying on the 
boundaries by constructing educated guesses based on several 
special classes of states that are simple to parameterize. 

A. Pure states 

We expand the pure states of two qubits in their Schmidt 
basis and write them as \iff)AB - 2^=1 V^IOa I*)b. where {|!)_ 4 ) 
and IIOb) are orthonormal basis of A and B and ci H- C 2 = 1, so 
they are easily parameterized by a single number ci. Thanks 
to this, we can easily show that 

I^<\0ab) - ^l‘^<\>^'}AB) = ^(1 - 2ci)2. (47) 

Moreover, Eq. (El 3) can be greatly simplified for pure states 
and the variance can be computed analytically, 

AJ---(I<A>ab)= ^(l-2ci)^ (48) 

In the end, we find the simple relation 

VAJ-<|(^>ab) = -^{r^mAB) - m^AB))- (49) 

These states provide the lower boundary for the right plot of 
Eig. 2 and are highlighted with a red dot-dashed line. A red 
dot-dashed line is also shown in the left plot for comparison. 
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Figure 2. (Color online) Left: Color plot of the square root of the variance of the skew information for 10^ randomly generated two-qubit 
states (gray dots), as a function of the AvSk and of the LQU. Right: Square root of the variance of the skew information plotted as a function 
of AvSk — LQU. Special classes of states (detailed in the main text) are highlighted by different lines, using the same style in both figures for 
comparison. 


B. Product states 

Next we consider the product states, which all have zero 
LQU. We seek a family of product states depending on only 
one parameter and interpolating between one state of the form 
( 1 , 4 / 2 ) ® pb, which has zero AvSk, and one pure product state, 
which has the highest AvSk among product states. Therefore, 
we consider the family of states given by 

Pprod(f’) := [p |0)a <0| + (1 - ) ® |0)b <0|, with p 6 [0,1]. 

(50) 

Their AvSk is easily computed 

-^‘^^(ppi-od(p)) " ^ 

Eq. (El 3) can be again evaluated analytically in this case and 
becomes 

AP^*^(pp,od(p)) = ^ (1 - ■ (52) 

In the end, we find the same simple relation as for pure states, 

^AJ-.(pp„d(p)) = ^(i'‘"'(pp™d(p))-'^/"^(Pprod(p))). (53) 

Therefore, these product states also lie on the lower boundary 
of the right plot of Eig. 2 and are highlighted with a green dot- 
dot-dashed line. A green dot-dot-dashed line is also shown in 
the left plot for comparison. 

C. Pure-QC states 

Another interesting class of states is given by the (pure 
quantum)-classical states, introduced above (see Eq. (36) in 


Sec. V B). Eor two qubits, these are states of the form 

PpQC = pI'Ao>a<<AoI®|0)b<0|h-(1-/7)|^i)a<i/'iI®|1)b<1|, (54) 

with |(/^o)a ^nd |iAi)a arbitrary pure states. A special subset of 
pure-QC states is obtained by taking 

PpQc(p) = ^|0)a<0| ® |0>B<0| + ^|+)a<+I ® IDiidI, 

(55) 

with p 6 [0,1], which linearly interpolates between a pure 
product state (when p = 1) and the maximally discordant sep¬ 
arable state of Eq. (45), i.e. the one having LQU equal to 
1 /2 [9] (when p = 0). These states have constant AvSk, with 
value 2/3. Their LQU and their variance can also be explic¬ 
itly computed as functions of p, although we do not report 
here the expressions. Combining the three quantities, we get 
the following relation 

^AJ-^(ppqc(p)) (56) 

1 + 3 J(j-*-(ppQc(p)) - 'l/-=(PpQc(p))) - . 

We conjecture that pure-QC states provide the upper boundary 
for the right plot of Eig. 2 as highlighted with a black dashed 
line. This is well supported by the numerical evidence. A 
black dashed line is also shown in the left plot for comparison. 

D. Separable states 

Erom Sec. VIIC we see that the state which behaves most 
differently (in terms of the variance) with respect to pure and 
product states is given by the maximally discordant separable 
state of Eq. (45). The leftmost upper curve for the right plot 
of Eig. 2 (shown by a black dotted line) connects this state p 
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with a state of the form (Ia/2) ® Pb- We then make an ansatz 
that separable states of the form 

p.sep(p) = I (|0>A<0| ® |0>B<0| + |+)a<+| ® |1>b<1|) + (X-p) 

(57) 

with p 6 [0,1], will attain the boundary. We see that once 
again our ansatz is well supported by the numerics. A black 
dotted line is also shown in the left plot for comparison. For 
states of the form given in Eq. (57), the relation between the 
variance, the average and the LQU is given by 

^AJ-.(p,ep(p)) = ^(J‘"^(psep(p)) - '7/‘"‘-(Psep(p))). (58) 


VIII. THE ROLE OF CORRELATIONS 


In this last Section we are going to discuss the influence 
played by correlations on the average and the variance of the 
skew information. 


A. Bounds on quantum correlations 


We have already seen that the amount of quantum (discord- 
type) correlations has no specific effect on the AvSk. How¬ 
ever, since the minimum susceptibility (i.e. the LQU) is in¬ 
stead a proper measure of quantum correlations, the variance 
of the skew information is affected as well. Roughly speaking, 
we can see that if two states have the same AvSk but one state 
has more discord-type correlations (i.e. higher LQU), its vari¬ 
ance will be smaller compared to the other. We can turn this 
into a quantitative statement and derive bounds for the LQU 
by combining the information about the average and the vari¬ 
ance of the skew information. In the case of two qubits, these 
bounds read 


max{0,/"^^(p) - ^|5M°■Xp)] < ^l‘"Xp) 

^ |j-=(p) - \ yjS M-Xp) 

- |j-=(p) - 7 - ^15AJ-.-(p) - i 


(59) 

if A J^-'^lp) <1/45; 
if 1/45 < AJ°-=(p). 


We expect that these bounds become tighter and tighter by in¬ 
troducing higher moments of the statistics. Indeed, knowing 
more about the distribution should also give more information 
about the minimum. For two-qubit states this is not very use¬ 
ful, since we can easily compute the LQU directly [6]. How¬ 
ever, in higher dimensions it is not known how to perform 
analytically the optimization involved in the computation of 
the LQU and one has to resort to numerical analysis. The 
approach presented here exploits quantities (the average and 
the variance of the skew information) that are exactly com¬ 
putable in any dimension (although their expressions could be 
rather involved) and could be then easily generalized beyond 
the two-qubit case. 


B. Connections with a recent measure introduced by Luo et al. 


We now show that the AvSk can be corrected to yield a 
proper measure of correlations. In order to do this we recall 
that a big contribution to the value of the AvSk comes from 
the purity of the local state of subsystem A. Roughly speak¬ 
ing, this accounts for the fact that a pure state \iIja) is more 
efficient than a mixed state (e.g. the maximally mixed state 
Pa - IaINa) in detecting the action of local operations. More¬ 
over, if we think in terms of the discrimination protocol de¬ 
scribed in [9], we can consider the situation in which we use 
a bipartite probe but perform only local measurements on A 
to extract the information (e.g. if we lose access to the refer¬ 
ence subsystem B). In the latter case, the only resource that 
we can exploit is the local purity of the state pA - '^'CbIpab}- 
Every benefit that we gain by measuring the global state pab 
must hence come from correlations. Motivated by the above 
reasoning, we define the quantity 


itXr^PAB) = I'^^ipAB) - I^^ipA ® Pb) 


Tr[A^] 

Nl-1 


(TrA[VTrB[PAB]]f 


(60) 

TrB[(TrA[ VpabI)^] ■ 


Note that the quantity I^^(pA ® Pb) in Eq. (60) actually de¬ 
pends only on pa and not on the other marginal ps = TrA[pAB]- 
We now show that the above quantity is equal (up to a pref¬ 
actor) to the quantity introduced by Luo and collaborators 
in [21]. They define 


F(pab) = Qa(pab) - Qa(pa 'S>Pb), (61) 


where Qa(p) is the average of the skew information with re¬ 
spect to any family of orthonormal hermitian operators, 
i.e. with respect to any orthonormal basis for the real Hilbert 
space L{9Xa) according to the scalar product (A, B) = Tr[AB]. 
That is, 


Qa(p)^Yj ^ Tr[A,A,] = d,,-. 

/=1 

(62) 

The quantity Qa(p) can be evaluated by writing the skew in- 

formation as in Eq. (4) and noting that A",® A,- = 5 a|a' IN a, 
where S a\a' is the swap operator. After some manipulations, 
the final expression Qa(pab) = Na - TrB[(TrA[ a/PabJ)^] can 

be found. We see that Qa{pab) coincides with I^^(pab) apart 
from a numerical prefactor that depends only on the choice of 
the spectrum. Therefore I^Xrr(PAB) is proportional to F(pab) 
in general. Luo et al. have shown [21] that F(pab) satisfies 
the following properties; 

- F(pab) = 0 if and only if pab ^ Pa® Pb- 

- F(pab) = F(Ua ® VaPabU\ ® yj) is invariant under local 
unitary operations, V Ua, Vb- 

- F{pab) is decreasing under arbitrary CPTP maps on sub¬ 
system B. It is also conjectured that F(pab) is decreasing un¬ 
der arbitrary CPTP maps on subsystem A. 

From these, they argue that F{pab) is a measure of total 
correlations, but cannot be specifically regarded as a measure 
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of classical or quantum correlations. We conclude that the 
same holds for our quantity I^orriPAB) defined in Eq. (60). 

Our analysis, though, complements the results of Luo et 
al. in two key points. First, we have provided a closed and 
compact expression that can be evaluated for any dimension of 
the Hilbert spaces 'Ha and 'Hb and that is not explicitly found 
in [21]. Second, we have given a clear operative interpretation 
to I^^„(pab), as the advantage that we can gain in doing some 
metrology task (e.g. state discrimination) by fully exploiting 
the amount of correlations in a bipartite probe. 

IX. SKEW INFORMATION AND METROLOGY 


we fix a set of isospectral generators //(A) for the unitary 
phase transformation the LQU [6] and the AvSk give 

strong indications about the minimum and the average estima¬ 
tion precision with respect to this set. 

Second, it was shown in [9] that the quantity 1 - Q(p,p') 
shares strong connections with the skew information. Specif¬ 
ically I(p, H) can be seen as the efficiency of a discrimination 
process where the two states that need to be distinguished are 
given by [p, e^'^pe‘^} and the unitary transformation is a 
small perturbation of the identity operator. That is, we have 
the relation 

lip, H )1 - Qip, e-‘^pe‘^) (64) 


As we discussed in the introduction and mentioned several 
times throughout the paper, the choice of the skew informa¬ 
tion has the twofold advantage of allowing easy manipulations 
and retaining interesting connections with the field of quan¬ 
tum metrology. Here we make these connections explicit, in 
the hope of conveying a clearer message to the reader before 
moving to the conclusions. 

Two common problems that are studied in quantum metrol¬ 
ogy are phase estimation [1-5] and state discrimination [1, 
36, 37]. In phase estimation, the goal is estimating a con- 
tinuos parameter (p that characterises the unitary transforma¬ 
tion p e~‘^‘^pe'^‘^ of the initial state of the probe. After 
choosing the best possible measurement strategy and the best 
possible estimator ptest (i-C- a function of the measurement 
outcomes and probabilities that is used to guess the actual 
value (p), the achievable precision in the limit of n » 1 repe¬ 
titions of the protocol is determined by the quantum Cramer- 
Rao bound [1, 38], which relates the variance of the estima¬ 
tor A^pbest "= il(nF(p, H)) to the inverse of the quantum 
Fisher information F(p, H). In state discrimination, the goal 
is discriminating between the initial state of the probe p and a 
transformed state p'. Since p and p' are in general not orthogo¬ 
nal, therefore not perfectly distinguishable, the figure of merit 
used in this case is the probability p of guessing correctly, 
that scales asymptotically in the number n of repetitions as 
p(n) 1 - Q(p,p'T, where Q(p,p') - min^ Tr[p^p''~^] is 
the quantum Chernolf bound [36, 37]. 

The skew information is strictly related to both the quan¬ 
tum Fisher information and to the quantum Chernolf bound, 
as we show in the following. First of all, the original quantum 
Fisher information [1, 38] is not the only generalization of the 
classical Fisher information, but there is a whole family of 
so-called generalized quantum Fisher informations [19, 20]. 
They all share a set of fundamental properties, e.g. are con¬ 
vex and have the same value for pure states, and they are all 
upper bounded by the original quantum Fisher information. 
The Wigner-Yanase skew information (multiplied by a factor 
4) belongs to this family [20] and this fact combined with an¬ 
other recent result [39], allows us to write 


1 


M\p,H] 


A 2 ”7 

— (pbest ~ 


1 


1 


F\p,H] - H\p,H] 


(63) 


We see that the skew information can be used to set upper 
and lower bounds to the estimation precision. Therefore, if 


where [/i, ] are the eigenvalues of H (see [9] for a formal char¬ 
acterisation). Moreover, the relation between the skew infor¬ 
mation and the quantum Chernoff bound becomes even more 
stringent when H is any operator acting on the Hilbert space 
of a qubit: indeed, in this special case the two quantities are 
proportional and we get 

lip, H) 1 _ g(p, e-'"pe'«). (65) 

We see that the LQU and the AvSk can therefore be used to 
characterise the minimum and average efficiency in discrimi¬ 
nating the elements of any of the couples |p, e^‘^'pe'^' ], ..., 
[p, where H\,... ,Hn belong to a set of isospectral 

Hamiltonians. 

All the above discussion remains valid even if we assume 
that p - pab is a bipartite state and the transformations 
act only on subsystem A, as we did throughout the paper. 
Moreover, with these additional assumptions we can use the 
LQU and the AvSk to draw another bridge between quan¬ 
tum metrology and quantum information theory, analysing the 
role of several resources in enhancing the metrological perfor¬ 
mance of different quantum states. It’s precisely in this sense 
that the analysis of Sec. VI acquires a strong relevance with 
respect to metrological applications. 


X. CONCLUSIONS 

In this work we have investigated a question of joint funda¬ 
mental and practical relevance, namely which resources and 
which bipartite states are useful as versatile probes, to achieve 
a required average performance in quantum metrology tasks 
[3, 5] involving a variable set of operations encoding an un¬ 
known parameter on one subsystem. We demonstrated that 
the average susceptibility of a state to local unitary transfor¬ 
mations can be reliably quantified by the average skew infor¬ 
mation, a quantity that we introduce and calculate in closed 
form for bipartite quantum states of arbitrary dimension. The 
average skew information is found to be a convex measure 
strongly dependent on the local purity of the probing subsys¬ 
tem, and in general requiring entanglement to reach its max¬ 
imum value. However, separable or even product states can 
still achieve fairly satisfactory degrees of average skew in¬ 
formation, meaning that they can be reliable metrological re¬ 
sources on average, when entanglement is not available. 
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The results of our analysis have been contrasted with the 
related, but different setting, in which the worst-case (rather 
than the average) susceptibility to local transformations is 
studied for bipartite states [6-10]. Such a worst-case perfor¬ 
mance can be quantified by the minimum skew information, 
known as local quantum uncertainty [6], which is instead de¬ 
termined entirely by quantum correlations of the discord type. 
By analysing comparatively the minimum, the average, and 
the variance of the skew information, we have identified the 
role of state purity, separability, and correlations to identify 
probe states with extremal properties, classifying their broad 
potential for metrological tasks such as parameter estimation 
and state discrimination. The general analysis has been illus¬ 
trated in particular in the simplest yet particularly relevant in¬ 
stance of two-qubit probe states, for which we have provided 
a complete numerical characterization. 

In this paper we were not concerned with another impor¬ 
tant issue in quantum metrology, i.e. how the precision of the 
estimation scales as we increase the number of “constituents” 
in the probe. Much is already known on the problem. For ex¬ 
ample, as mentioned in the introduction, one can show that by 
using pure entangled states of n qubits the minimum estima¬ 
tion error can be reduced by a factor ^[n with respect to using 
a pure separable state [3, 5]. Moreover, a recent work [14] has 
provided evidence that a similar enhancement can be found 
for discordant mixed states of n qubits over classical mixed 
states, under particular measurement strategies. We remark 
that in this paper we didn’t study the role of correlations (and 
other properties of the probe) with respect to optimal perfor¬ 
mances. Instead we focused on a complementary aspect, i.e. 
versatility, and left outside, at least for the moment, consider¬ 
ations regarding the scaling of our functionals. 

We expect that the study of the minimum and average skew 
information for continuous variable systems [40-42] would 
provide us with further insights In this case, for example, 
it comes naturally that one does not have experimental ac¬ 
cess to the whole infinite-dimensional Hilbert space and must 
work with limited resources (e.g. limited classes of states 
and operations, limited energy, limited squeezing, limited 
amount of entanglement, limited purity, ...). The minimum 
and the average skew information would provide then clear 
and simple-to-evaluate criteria that can help in picking opti¬ 
mal probes among the set of accessible states. Based on the 
recent progress in calculating some of these measures (such 
as the interferometric power and the discriminating strength) 
for Gaussian states of continuous variable systems in a worst- 
case scenario [8, 10, 43], we believe that a Gaussian version 
of the average skew information might be amenable to analyt¬ 
ical evaluation; it would then become particularly important 
to study its scaling with the resources typically involved in 
optical interferometry, such as the mean energy of the prob¬ 
ing system [8], and with other nonclassical features such as 
squeezing and entanglement. This is left for future investiga¬ 
tion. 

This work has provided yet another application of the 
Wigner-Yanase skew information, defined more than half a 
century ago [17], in quantum information theory. The skew 
information represents one of the most insightful and mathe¬ 


matically convenient quantum generalizations of the classical 
Fisher information [19,20], and it has proven useful already to 
derive improved uncertainty relations [18], to define measures 
of asymmetry (coherence) [44, 45] and correlations [6, 21], 
and to construct generalized geometric quantum speed lim¬ 
its [46]. The latter application, in particular, deals with the 
question: How fast can a quantum state evolve under a closed 
or open system evolution? The study presented in this work 
can be framed in a similar perspective, as the average skew 
information introduced here quantifies precisely how fast, on 
average, a quantum state of a bipartite system evolves under 
any local unitary dynamics (within a fixed spectral class) af¬ 
fecting one of its subsystems. The more versatile probes for 
quantum metrology are exactly those whose reaction to the 
local dynamics is faster, indicating an increased susceptibility 
to the unknown parameter encoded in the dynamics itself 

It is finally interesting to comment on the information- 
theoretic resource unlocking such an enhanced susceptibil¬ 
ity to local dynamics. If the figure of merit is the minimum 
susceptibility, the resource is local asymmetry (coherence) in 
all possible reference bases for the probing subsystem, which 
is equivalent to discord-type quantum correlations [6, 7]. If 
the figure of merit is the average susceptibility, instead, we 
demonstrated that the resource is local purity for the probing 
subsystem. This is clear when one considers that the “free” 
states with vanishing average skew information are those of 
Eq. (8), taking the form of a product of the maximally mixed 
state for the probing subsystem, tensor any state for the other 
reference subsystem. Therefore any degree of local purity be¬ 
comes useful in this context. This suggests that the average 
skew information could be further investigated as a quantum 
thermodynamical resource [47]. Namely, considering the case 
in which the probing system A has all degenerate energy lev¬ 
els (so that the maximally mixed local states are the only free 
states), the average skew information defined in this paper 
might be related to the amount of work that can be extracted 
from A by some optimal thermal machine with access to the 
reference storage system B, provided the machine is coupled 
to a heat bath [48]. Investigating these intriguing connections 
further will be the subject of an independent study. 
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Appendix A: The swap operator 

Let 'Hx and ^x' be isomorphic Hilbert spaces spanned by the orthonormal basis {|!)x) and {|i)z')- The swap or flip operator 
is defined by the relation 5x|x'|i)xlj)x' = \j)x\i)x' [26]. A possible representation is given by 



(Al) 


We report some useful properties that we use throughout the paper. 

1- (Sx\X')^ - txs where \xx' = Ix ® Ix' is the identity operator on “Kx ® 'Hx '; 

2. Sx\x’ = 5a|a' if^x ='Kr and'Kx' ='Ka' ® “Kb', with'TY a isomorphic to'K a' and'K b isomorphic to'K b-; 

3. 5x|x'(0x ® ^^x')‘5x|x' = (f^x ® ©x')j where 0x,x' and Ox,x' are linear operators on the corresponding Hilbert spaces; 

4. (0x ® Qx')‘5x|X' = ‘5x|x'(f^x ® 0X'). which simply follows from Property 3 if we apply 5x|X' to both terms; 


5. Trx,x'[(0x ® i^x')‘5x|x'] - Trx[0x^^x]- 


For completeness, we sketch the proof of Property 5. Without loss of generality we set 



(A2) 


By explicit computations we have 


Trx[0x^^x] - > 


(A3) 


ij 


and 


Trx,x' [(0x ® ^X' )S x|X' 1 Oija)emTv[\i)x{j\a)x(fi\ ® K)x'<'«lyS>x'<a'l] 


ij t,m 0',y8 



(A4) 


thus concluding the proof. 

Appendix B: The twirling channel 

Let 'Kx and ^x' be isomorphic Hilbert spaces and 0xx' an operator acting on the tensor of the two 'Kx ® ^x' ■ The twirling 
channel modifies this operator by applying the same local unitary operation simultaneously to X and X' and then averaging this 
action over all possible local unitaries. This is expressed as 



(Bl) 


where dpniU) is the Haar measure over the unitary group {U{Nx)] and Nx - dim('Kx). It can be shown that the integral in 
Eq. (B 1) has a simple solution in terms of the swap operator as introduced in Appendix A 



(B2) 


Appendix C: Isotropic states of two qubits 


We consider the isotropic states of two qubits [27] 


l-F 


AF - 1 


Pf - 


3 


3 


l<A+)(<A+l> 


(Cl) 
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parametrized by 0 < f < 1, where F is the fidelity between the isotropic state and the Bell state \ij/+) - (|00) + |11))/ V2. The 
AvSk can be easily computed by decomposing the identity on the Bell basis. 


Pf 


l-F 

3 




l-F. 




l-F 

3 


+ Fli/r+Xl/r+l, 


(C2) 


SO that 


vTt - 


l-F. 




1-F. 


-|0+><<^+ 


l-F 




(C3) 


Plugging this into Eq. (20), we get that the AvSk is given by 1 


. Similarly, the LQU can be computed 


following the prescription of [6] and the result is again given by the same expression 1 - 


, as anticipated 


in the main text. Isotropic states of two qubits lie on the blue line in Fig. 1 . Entangled isotropic states (with F > 1/2) have 
AvSk=LQU> ~ 0.09. Therefore, this gives a simple proof that some entangled states lie on the left of the curve LQU= 1 /2. 


Appendix D: Integrals over the unitary group 


Consider a general integral of the form 


f 


^hji • 




(Dl) 


where t/ is a unitary matrix acting on the Hilbert space Fix, dfiniU) is the Haar measure over the unitary group {U(Nx)} and 
Nx - dim('?/x). Such integral is called a moment of order (n, m) of the unitary group and is zero whenever n F m [49]. In the 
case n - m, the integral can be computed using the Weingarten calculus [50] and yields the expression [5 1 ] 


n 


n ^r(n) ^ ja ^Ta-{a) ’ 

(TjeSn a=\ 


(D2) 


where cr and t belong to the symmetric group, i.e. they are permutations of n elements, and c(n, cr) are the so called Weingarten 
functions which depend on the number of elements appearing in the integral and on the particular permutation of those n 
elements. The analytic expression of the Weingarten functions is explicitely known for small values of n [51], and it can be 
computed for higher n with some effort. 

Note that the twirling channel implicitly contains an integral of the form (Dl), with n - m-2. 


1. Weingarten functions for the case n = m = 4 


The case n = m = 4 is particularly interesting to us since it appears in the computation of the variance of the skew information. 
We report here the Weingarten functions for n - 4, taking them from [51]. First let us set some notation to deal with permutations. 
A permutation cr of 4 elements will be written as its action on the string (1,2,3,4]. So for example, the permutation (3 2 4 1) 
maps (1,2,3,4] to (3,2,4,1], i.e. brings the first element to the fourth place, the third element to the first place, the fourth 
element to the third place and leaves the second element unchanged. Where possible we will index these permutation by the 
associated permutation class, e.g. (3 24 1)^ [cr] = [1,3] as given by one cycle over 3 elements (1,3 and 4 in the example) and 
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one cycle over one element (2 in the example). With this in mind, we can now write 


cr 

c(4, cr, Nx) 

[cr] - [4] 

-5 

{Nx - 3){Nx - 2){Nx - l)Nx{Nx + l)(A^x + 2)(Nx + 3) 

M = [l,3] 

2NI-3 

{Nx - 3){Nx - 2){Nx - l)Nx{Nx + l){Nx + 2){Nx + 3) 

[cr] = [22] 

N^ +6 

{Nx - 3){Nx - 2 ){Nx - l)Nx{Nx + l){Nx + 2){Nx + 3) 

[cr] = [l2,2] 

1 

{Nx - 3){Nx - l)Nx{Nx + l){Nx + 3) 

M = [11 

N^ - 8^2 H- 6 

{Nx - 3){Nx - 2 ){Nx - l)Nx{Nx + l){Nx + 2){Nx + 3) 


In other words, the Weingarten functions depend only on the class of the permutation. An important thing to notice is that the 
Weingarten functions with n - A diverge for Nx < 4 (and in general they diverge when Nx < n). However, it has been proven 
that the sum in Eq. (D2) does not diverge because the poles in each term cancel out after careful simplifications [49]. This allows 
us to compute the variance of the skew information even for two-qubit states (having Nx = 2). 


Appendix E: Second moment of the skew information for two-qubit states 


The second moment of the skew information is defined as 




(El) 


This expression can be expanded as 


Tr[ xI^UaAIu{ + Tr[ x/^UaNaUI x/^UANAUlf 

- 2 Tr[ yf^UAKAUlUA^AUl ^/^] Tr[ ^UaAaU\]). 


(E2) 


The first term in Eq. (E2) can be integrated using the properties of the twirling channel and the result reads 


I 


2,,t.^n2_2Tr[A^]^-Tr[A-] ^ 2Tr[Al] - TilA^f^^ 


cI^hCUa) Tr[ ^[pUaA^UI yfp] 


-Tr,4[(TrBM)1. 


(E3) 


The second term in Eq. (E2) can be rewritten by expanding each operator on a basis {|/),4 \j}g} of 'TTas- For simplicity we fix 
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the local basis of A to make diagonal. 


df^HiUA) Tr[ ^JpUAAAU\ ylpUA^AU^f = 
J" dflH(UA)[ ^ ^/p hji,‘ 


hji ^hh-^h ^ UAa ^ u 


X y yfpi' /'/■' Ui- i' A/' Ul yfpi- i' i' r Ui' i'Ki' Ul 

/ i yr'\Jv‘‘2J2 '2’'3 '3 ’^'4-/2’'5’-/l '5’'6 'h 

i'f 

y VPhjiAji ^UhM,ii 

X r duH{UA)(Ui, i, Uu L Ui’ i' Ui' i' u^. . t/:''. ul., ul ). 

I r-n\ 12,12 15,<6 i2'h ^3,4 4,4 /g,/!/ 




Using the results of Appendix D, we can compute the integral over the unitary matrices and Eq. (E4) reduces to 


y c(4, cr, 2) A(Tcr, A) G 2 (t,p) 


(E4) 


(E5) 


where 


A(to-,Aa)= y AijAj^A/'Aj-^i 

44^3^6 


l3,l6,'3.'6).'rAI'3,i6,i3,'6)) 


(E6) 


is a function which depends only on the spectrum and on the composition of permutations tct, while 




hhU's iWVs jdi j\f2 


(E7) 


is a function which depends only on the state and on the permutation r. Both F and G 2 can be analytically computed for each 
choice of cr and t. Assuming that the spectrum is traceless, we find that F depends only on the class on the permutation tct 


TCT 

F{t(t, Ka) 

[tct] = [4] 

Tr[A4] 

[Ttr] = [1,3] 

0 

[rtr] = [22] 

Tr[A2]2 

[tct] = [12,2] 

0 

[tct] = [14] 

0 


while G 2 assumes different values even among permutations belonging to the same class 


T 

G2(t,p) 

T 

G2(t,p) 

r 

G2(t,p) 

(1 2 3 4) 

1 

(2 3 14) 

£ 

(3 4 12) 

D 

(1 2 4 3) 


(2 3 4 1) 

T 

(3 4 2 1) 

@ 

(1 3 2 4) 


(2 4 13) 

T 

(4 12 3) 

r 

(1 3 4 2) 

£ 

(2 4 3 1) 

£ 

(4 13 2) 

£ 

(1 4 2 3) 

£ 

(3 12 4) 

£ 

(4 2 13) 

£ 

(1 4 3 2) 


(3 14 2) 

T 

(4 2 3 1) 

S 

(2 13 4) 


(3 2 14) 

S 

(4 3 12) 

@ 

(2 14 3) 

C 

(3 2 4 1) 

£ 

(4 3 2 1) 

D 


(E8) 


(E9) 
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The various letters in the last table are shorthand notation for the following expressions 

S = Tr^ [(TrB[p])2], 

D = Tr^A' [(TrB[(lA- ® -y^) • (Ia ® yfpA^)]f\ , 

£ = Tr[v^■(TrB[p]®TrA[^4^])], 

T = TrA [(TrB[ Vp ■ (Ia ® TrA[ , 

Q - Ti'aa' [(TrB[(lA' ® ^Ipab) ■ (Ia ® VP^'b)])^‘5a|A'] • (ElO) 


In the last expression, S a\A' is the swap operator [26] discussed in Appendix A. 

The third term in Eq. (E2) can be tackled similarly to the second term (note that we explicitly wrote it with Ua and t/^ 
appearing 4 times each). The result is 

-2 Yj c(4, 0 -,2) FCtct, A) G3(t,p) (Ell) 


where F is the defined in Eq. (E 8 ) and G 3 is defined below 


T 

G3(t,p) 

T 

G3(t,p) 

T 

G3(t,p) 

(1 2 3 4) 

1 

(2 3 14) 

1 

(3 4 12) 

£ 

(1 2 4 3) 


(2 3 4 1) 

Jl 

(3 4 2 1) 

S 

(1 3 2 4) 


(2 4 13) 

Jl 

(4 12 3) 

2£ 

(1 3 4 2) 

£ 

(2 4 3 1) 

Jl 

(4 13 2) 

2B 

(1 4 2 3) 

£ 

(3 12 4) 

2S 

(4 2 13) 

£ 

(1 4 3 2) 

S 

(3 14 2) 

2£ 

(4 2 3 1) 

S 

(2 13 4) 

2 

(3 2 14) 

S 

(4 3 12) 

S 

(2 14 3) 


(3 2 4 1) 

£ 

(4 3 2 1) 

£ 


By putting together Eqs. (E3), (E5) and (El 1) we find 

, . . _ 2Tr[Aif-Tr[Al] 2Tr[Al] - Tr[Ai]^ 

\ '^V(£/p 6 6 

+ 2 c(4, cr, 2) F{Tcr, A) (G 2 (t,p) - 2 G 3 (t,p)). (El3) 


Unfortunately, we cannot further simplify this expression to explicitly show that the poles appearing in each c(4, cr, 2) cancel 
out. However, a direct computation proves that Eq. (El 3) does not diverge. 
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